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Abstract 

The ground state of a many body Hamiltonian considered in the quasiparticle representation 
is redefined by accounting for the quasiparticle quadrupole pairing interaction. The residual in¬ 
teraction of the newly defined quasiparticles is treated by the QRPA. Solutions of the resulting 
equations exhibit specific features. In particular, there is no interaction strength where the first 
root is vanishing. A comparison with other renormalization methods is presented. 
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I. INTRODUCTION 


fl 


The big merit of the liquid drop model (LDM) proposed by Bohr and Mottelson jjJ is that 
one defined by the concept of rotational bands. Also, some collective properties of spherical 
nuclei have been nicely described. The main drawback of LDM consists of the fact that it 
accounts only for the spherical and harmonic motion of the drop, while many experimental 
data reclaims a non-harmonic picture and, moreover, many nuclei exhibit static deformed 
shapes. Many phenomenological improvements have been proposed along the time, among 
which few are to be mentioned: a) rotation-vibration model jsj]; b) Gneus-Greiner model Q ; 


approximation 


Q. 


c) generalized collective formalism 0 ; d) coherent state model |5|, 6]; e) interacting boson 


n 


In parallel, microscopic theories have been formulated trying to get 


counterparts of the phenomenological methods and interpret the nuclear collective motion 
in terms of the single particle motion. Thus, the random phase approximation built on the 
top of either the Hartree-Fock or the BCS ground state (QRPA) [gj] provides a collective state 
which corresponds actually to the one phonon state predicted by the harmonic LDM |q]. 
Another important result is that of Kumar and Baranger, who calculated the inertial and 
stiffness parameter microscopically jl^ t the potential energy surface leading to some sound 
nuclear structure interpretation. Based on the RPA ground state, several procedures of 
accounting for some new correlations, i.e. of going beyond RPA, have been proposed. Such 
procedures are related with the equations of motion method 
technique Q H, Q, Q. 

A distinct procedure which keeps the appealing harmonic picture of RPA but also includes 
in the definition of the phonon operator new correlations is obtained by renormalizing the 
specific equations 20]. Thus, the drawback of the standard RPA formalism of collapsing 
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14] or boson expansion 


for a critical value of the attractive long range interaction strength is removed. Indeed, the 
collective root of the RPA equations goes to zero not for a finite value of the mentioned 
interaction strength, but only asymptotically. This approach was extended to the proton- 
neutron dipole interaction in Refs. j^lj]. Vanishing the excitation energy of the collective 
RPA state corresponds to a phase transition where the ground state is unstable to adding 
small contribution. Around this critical interaction the RPA method is no longer valid. In 
order to stabilize the ground state it is necessary either to change the mean field for the single 
particle motion, which results in having deformed single particle orbits, or to renormalize 
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the basic equations. 

In the present paper we propose a new method of renormalizing the QRPA equations. 
As we shall see, the result for the collective root is that it does not vanish in the critical 
interaction strength, where the standard QRPA collapses, but reaches a minimum value and 
moreover the energy increases when a subsequent increase of the strength is performed. The 
new point of this work is that the mean held is redefined in the quasiparticle picture by 
including in the ground state the quasiparticle quadrupolc pairing correlations. As a result, 
both the new quasiparticles and the new QRPA solutions are deformed. 

The project sketched above will be described according to the following plan. In Section 
2 the model Hamiltonian is presented. For the sake of completeness the results for the 
standard BCS and QRPA equations are briefly described. Section 3 is devoted to the 
deformed quasiparticles or, in other words, to the second order BCS approach. The new 
BCS and QRPA equations are analytically derived. In section 4 the formalism is numerically 
applied to a single j shell. The final conclusions are drawn in Section 5. 

II. THE MODEL HAMILTONIAN 

We consider a system of nucleons described by a many body Hamiltonian consisting of a 
spherical mean held term, the pairing and the qnadrnpolc-qnadrnpolc two body interactions. 
Written in second quantization, this has the form: 

= 5>« - A )4e„ - jP'P -jY. QvQi-A-r, (2.1) 

a fi 

where Pt and Q denote the pairing and qnadrnpole operator, respectively: 

p' = 

a 

= ^(a|r 2 F 2M |/3)4c/3 = ^ q ab (4q>) 2m , 

a,(3 a,b,fi 

Qab = i(a||r 2 y 2M ||Q, with j a = \/2ja + 1. (2.2) 

The particle-hole qnadrnpole operator is defined as: 

(4q,) 2 „ = E ■ (2.3) 

m a ,mp 
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The second quantization is used within the spherical shell model basis |a) = | a,m a ) = 

I nJajaiTi Q ). Thus, the creation (annihilation) operator of one particle in the state la) is 
denoted by c* (c Q ). We used also the notation | — a) — |a, —m a ). 

The sum of the first two terms is quasi-diagonalized by passing to the quasiparticle 
representation defined by the Bogoliubov-Valatin (BY) transformation: 


fit c) v l / a'S Q ,C Q! , 

a a = U a c a - V a s a ct a , with s a = (-) J 


a rrra: 


(2.4) 


The quasiparticle vacuum state will be hereafter denoted by | BCS). Up to an additive 
constant, in the new representation the Hamiltonian is: 

H = ^Kaia„~^Y.Q^.-^-y< ( 2 . 5 ) 

a ij, 

where E a denotes the quasiparticle energy for the state characterized by the set of quantum 
numbers ”a”, while the quadrupole operators can be expressed in terms of two quasiparticle 
and quasiparticle quadrupole density operators: 


Q 21 J, 

A^ ah ) 

Bl(ab) 


Y Qo-b tab (Au,( ab ) + A 2 + Vab (^(06) + H 2 _ (U (o6)(-) 


a<b 


vT+i; 


ab 


E nd a j\ 

mocTT 


i b 2 n f r/ t 

mpp Li a u l3i 


,mp 




,mp 


A 2 p(ab) = (^(aft))*; B 2li {ab) = (B^iab)^, 

U = 7 ====== (U a V b + U b V a ), Vab = -A— (U a U b - V a V b ). (2.6) 

V 1 + 0 ab 1 + Oab 

In deriving the expression of the quadrupole operator in the quasiparticle representation, 
some symmetry properties were used: 


An(ba) = H Ja 3b Ap,(ab)] A 2fl (ba) = (-) 3a 3b A 2tl (ab), 

Bl(ba) = {-y^B 2 _,{ab)(-Y- B 2 ,(ba) = (-)*-*Bj _ #1 (a6)(-)'*, 

q ba = (- ) ja ~ jb q a b■ (2.7) 


The quasiparticle many body Hamiltonian is treated within the Random Phase Approx¬ 
imation (QRPA) formalism. Thus, one defines a phonon operator 


Cl = Y [x(ab)Al(ab) - Y (ab) A 2 ^(ab)(~y 


( 2 . 8 ) 
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with the amplitudes X(ab ) and Y(ab ) determined such that the following equations are 
fulfilled: 


H,C{ 


2/x 




c 2 „, cl, 


= S, 




(2.9) 


The first equation yields for the phonon amplitudes the so called QRPA equations: 


(a b 

\-B* -A* 



( 2 . 10 ) 


This is a homogeneous system of linear equations which determines the amplitudes up to a 
multiplicative factor which is fixed by the second equation (12.9(1 . which gives: 


53 [|V(a6)| 2 - |Y(af>)| 2 ] = 1. (2.11) 

ab 

The matrices involved in Eq. (12.9ft have the expressions: 

X 

Aab^a'b' ( E a T E^j b a a f Sfrb' ^ I'ah^a' 1/ • 

B a b,a'b' = -y r ab r a ' b ', with r ab = q ab £ ab . (2.12) 

Once the QRPA equations are solved, the phonon space is defined. Thus, the vacuum state 
denoted by | RPA) is the ground state, while the excited states are multi-phonon excita¬ 
tions of |RPA). Since the two body interaction is of a separable form, the QRPA equations 
provide a dispersion equation for the excitation energies and analytical expressions for the 
phonon amplitudes. By inspection of these expressions, one finds that the first excited state 
corresponds to an energy smaller than the minimal two quasiparticle energy and has a col¬ 
lective character. Due to the attractive nature of the two body interaction, the collective 
state energy is decreasing when the interaction strength X is increased. Consequently, for a 
critical value of X the phonon energy is vanishing and the phonon operator is undetermined. 
This situation defines the breaking down point of QRPA approach. To avoid this regime, the 
mean held for the single particle motion should be re-defined, which results in renormalizing 
the ground state energy such that the collective state has a positive energy. Concretely, 
the spherical shell model single particle basis is to be replaced with the Nilsson single par¬ 
ticle basis. Another way to remove the QRPA breaking down is to define a new phonon 
operator by going beyond the QRPA approach, which is actually based on the quasi-boson 
approximation of the two quasiparticle quadrupole operators A\ fl and A 2/i . The quasi-boson 
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commutation equations have been corrected by retaining from the exact expression not only 
the constant term, but also the scalar one which is considered in the average. This average is 
determined self-consistently with the QRPA solution and, consequently, the collective root 
energy goes to zero only asymptotically. In next section we formulate a new method to 
renormalize the ground state energy. 

III. DEFORMED QUASIPARTICLES 

Here we study the BCS ground state excitation: 

| BCS) = with 

T ZqCL^CL —q, ZqQj—qlCIolSol. (3.1) 

For what follows it is useful to employ the polar representation of the parameters z a \ 

z* = p a e itp . (3.2) 

The images of the quasiparticlc operators through the aforedefined transformation, are: 

dl = e T o) a e~ T = al cos(2 p a ) - a- a s a sm(2p a )e~ t<Pa , 
d a = e T a a e~ T = a„cos(2p Q ) — al a s Q sm(2p a )e lipa . (3.3) 

With the obvious notations: 

u a = cos( 2/9 q ); u a = sin(2 p a )e~ liPa , (3.4) 

we recognize the BV transformation for the quasiparticle operators. The transformation 
parameters u Q and v Q satisfy the equation: 

w a + |w| 2 = l, (3.5) 

which reflects the fermionic character of the new quasiparticle operators d\ and d a . We 
note that the transformed state | BCS) is vacuum state for the quasiparticle annihilation 
operators d a . It is worth mentioning that due to the transformation dependence on the 
magnetic quantum number rn a , the new quasiparticlc operators are tensors of rank j a with 
indefinite projection. In this respect, one can assert that the newly defined quasiparticles 
are deformed operators. 
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The parameters u a and v a may be viewed as classical coordinates depending on time. 
Moreover, considering the polar representation, the phase ip a , defining the coefficient v a , has 
the meaning of a generalized linear momentum. Consequently, since we are concerned with 
the static properties of the new BV transformation coefficients, it is reasonable to consider 
vanishing phases ip a . 

In what follows we try to determine the parameters u a and v a such that the transformed 
state | BCS) becomes the true ground state of the many body system under consideration, 
i.e. it corresponds to an energy lower than that associated with | BCS). The new ground 
state is a deformed function, which is reflected by the non-vanishing value of the expected 
quadrupole moment: 


{BCS\Q 20 \BCS) = V20 J2 '¥ctlL. + u a v. a ) +,Uvl + O) = <Zo- 

a,m a >0 ^ a 

(3.6) 

The last part of the above relation expresses the fact that the state | BCS) has a definite 
quadrupole deformation q 0 . The average value of H with the deformed state | BCS) is: 


S' = (BCS\H\BCS) = Y J E a vl- 


A 2 

^20 

AX 


(3.7) 


X 


r\ '/’ ^ \/la,b^'rn a —m a 0 (^o-b{ u —a u b,m ct v a v b,—m a )A-‘JJ a b{ u b,m a v a¥U— a V} ) — rna ))^ 


a,b,m c 


where we denoted: 

A 20 = ^-(BCS\Q 20 \BCS). (3.8) 

Note that the last term of Eq. (13.8p is provided by the average of the quasiparticle terms 
of the type d a db trna d\, m l d^_ a ,. However, such terms will be treated at the QRPA level, i.e. 
at a later stage. Due to this reason, hereafter, the mentioned term from Eq. (13. 8 p will be 
neglected. 

In what follows we shall look for the stationary points of the function: 

f = + 1). (3-9) 

a a 

with n a denoting the Lagrange multiplier corresponding to the restriction (13.51) . These are 
solutions of the equations obtained by vanishing the first derivatives of E with respect to 
the parameters u a and v a , defining the BV quasiparticle transformation. Eliminating the 
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Lagrange multipliers, one arrives at: 

A 2 0 CaaQaa (WqU—q WT—a) {E a 2A 2 o?7 aaQact) ^ aVa 0> 

where the following notation has been used: 

2 

f) — _ (~<ja2ja 

^Cota ~ a ^ m a 0m a Haa- 

2 Ja 


(3.10) 


(3.11) 


We remark that the above equation is invariant to the change a —> —a. This suggests that 
the solutions of the above equation satisfy: 


XLql OLi ^ot V— q . 


(3.12) 


Thus, the occupation probabilities acquire the expressions: 

E a — 2A 2 o?? aaQcx. 



l=i( 

\<l 

! 


sj(E a - 2A 20Vaa Q aa ) 2 + (2A 20 £ aaQaa ) 

It is interesting to note that by the quadrupole moment restriction and 

A = J°_ 

20 4A ’ 

Eq. (l3.13p fully determines the parameters v a and u a . With the notations: 

E a E a 2,/\ 2 oQ aon 2A 2 oQ aa , 


(3.13) 


(3.14) 


= \/El + 5l, 

one obtains a more transparent expression for the BY transformation coefficients: 


W i = 1 f 1 Ea 

u 2 I 2 \ e, 


(3.15) 


(3.16) 


The quantity e a has the significance of the second order quasiparticle energy, i.e. the energy 
corresponding to the deformed quasiparticle state |a) = d^BCS). Further, we shall define 
a phonon operator 


r f = 52 (*i 4 < 

ab,m a 

such that it obeys the equations: 


Y ab dt (I ^ 

1 m a u, o-m a u a I , 


(3.17) 


[ff,r*]=ur*, [r,r*]=i. 


(3.18) 
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In terms of the new quasiparticles, the model Hamiltonian is: 


H 




j E (44,+ d b ,- m j a ) 

a,b,m a 


E <?" (44,- m „, + ^,-^44) 

a',b',m a i 

(3.19) 


with the notation: 


Qm a — ^m ol -m a Q < iab [^(u a Ub,-m a ) + Vab ( u -aVb,-m a "P v aU b ,m a )] 


(3.20) 


The amplitudes X^ a and are determined by the QRPA equations, which are of a similar 
form as those given by Eq. (I2.10p . and the normalization condition: 


2 E [At) 2 - «) 2 

a,b,moc 


l. 


The matrices involved in the QRPA equations have the expressions: 


(3.21) 


A 

B 


ab;a'b' 

m a ]m a / 

ab;a'b' 

m a -m a , 


(e a + e b>ma )5aa'5 w 5 mama , - XQ^Q^, 

_ -y fxab r\a'b' 

A Vm a Vm a ,- 


(3.22) 


Since the two body interaction involved in H is separable, the compatibility condition for 
the QRPA equations may be brought to the form of a dispersion equation, while the phonon 
amplitudes are analytically expressed. Now, it is worth noting that both the QRPA and the 
deformed BCS equations involve the factors i] a b- This reflects the fact that the terms B 
and B 2 -^(—) m of the Hamiltonian expressed in terms of spherical quasiparticles, contribute 
to the mentioned equations. This feature contrasts the standard QRPA equations, which 
ignore the scattering terms . 


IV. THE CASE OF A SINGLE j SHELL 

The essential features of the QRPA formalism with a multi-shell calculations can be 
recovered by restricting the single particle space to a single j. Since here we are not interested 
in quantitative details, but rather in underlying the main virtues of the proposed formalism, 
we consider the numerical application for the single j case. Thus, one considers a system 
of N=10 nucleons moving in the spherical shell model state j = A 3/2 and described by 
the corresponding many body Hamiltonian (12.ip . We present separately the spherical and 
deformed QRPA results. 
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A. Results for QRPA built on the top of the first order BCS 


The occupation probabilities are: 


E 2 


N 

20 ’ 


U 2 



(4.1) 


For the sake of simplifying the notation, the low indices of U and V, specifying the chosen 
single j are omitted. The state semi-degeneracy is denoted by fl Neglecting the term 
GVtV 4 accounting for the renormalization of the single particle energy due to the residual 
interaction, the BCS ground state energy is: 

A 2 

E = 2eQV 2 - — = 2eQV 2 - GQ 2 V 2 {1 - V 2 ). (4.2) 

G 


Here the Fermi level energy is set equal to zero. The condition of minimum energy leads to: 


V 2 = - 1 


2e 
GQ 

This expression is consistent with Eq. (14.ip if the single particle energy is: 


(4.3) 


a si 


e = 


1 - 


N 

77 


(4.4) 


Equation (j4.3j) gives for the quasiparticle energy: 

GQ 


Eq = 
q 2 


(4.5) 


In our application we took G = 0.4MeH, which results of having E q = lAMeV. The 
compatibility condition for the QRPA equations reads: 


= 4q 2 - 2x44. 


(4.6) 


The positive root of this equation is: 


cj = 


G 2 d 2 ~\SIX ({jWr^Wj )) 2 A 


N \ 

l ~2n) 


1 V 2 


(4.7) 


The QRPA energy is plotted in Fig. 1 as function of X. From this figure we notice that 
for X = 0 the mode energy is equal to twice the quasiparticle energy, while for X & 
1.82 10~~ 3 MeVfm~ 4 this is vanishing. Within this interval the function is monotonically 
decreasing. The vanishing mode energy reclaims a breaking down for the QRPA approach. 
In next subsection we describe a method of recovering the validity of the QRPA approach. 
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10 5 X [ MeV fm 4 ] 


FIG. 1: The spherical QRPA energy as function of the quadrupole-quadrupole interaction strength 
for the case of a single shell, j = ii 3 / 2 - 

B. The study of QRPA for the quadrupole pairing correlated quasiparticles 

Switching on the quadrupole pairing correlations for the spherical quasiparticles, we have: 

E m = E,- j%(U 2 -V 2 )(j\\r 2 Y 2 \\j)C% m , 

S ,» = ^9cV2£/V'(i||r 2 y 2 ||i)Ci 2 L, (4.8) 

e <n — y O; + '() ■ 

vi = \(l-—), ui = l-vl (4.9) 

Numerical results correspond to q 0 = 20 fm 2 . One notices that the newly defined quasi- 
particle energy depends on the magnetic quantum number ”m”. The split, caused by the 
quadrupole moment of the single particle state | jm), is shown in Fig. 2 as function of X, 
the strength of the QQ interaction. Another peculiarity for the second order BCS ground 
state is that the average number of quasiparticles is not vanishing. 

(BCS\N,\BCS) = < 4 ' 10 > 


11 




FIG. 2: The deformed quasiparticle energies for the j = i 13 / 2 multiplet. 


where the quasiparticle number operator is denoted by 


N q = y(4.11) 

a 

The dependence of the quasiparticle average number on the interaction strength is shown in 
Fig. 3. Note that the larger the strength A", the larger the quasiparticle averaged number. 
This result implies also the presence of quasiparticles in the QRPA ground state. 

In order to write the QRPA equation, we need to know the matrix Q^ a defined by Eq. 
(13. 201) . In the case of a single j shell one obtains: 


Qii = (—r m c2 m (j||rV 2 ||j) (V2UV(u 2 m - v 2 J + (t/ 2 - V 2 )u m v„ 


(4.12) 


The compatibility condition for the QRPA equations can be written under the form of a 
dispersion equation: 




(OH)' 

4 p? — uj^ 

m> 0 771 


(4.13) 


For the chosen value of j there are seven solutions for u denoted by u(m) and ordered as: 


w(l) < u( 2) < ... < u( 7). 


(4.14) 


The first root has a collective character since several quasiparticlc pairs contribute to the 
phonon operator. Note that u;(l) is in magnitude smaller than the minimal two quasiparticle 
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10 4 X [MeV fm' 4 ] 


FIG. 3: The average number of quasiparticles in the second order BCS state, | BCS). 

energy, min(2e m ). The two aforecompared quantities are represented as function of the long 
range interaction strength X, in Fig.4. We note that cu(l) is no longer a monotonic function 
of X. There are, however, two intervals of different monotony. On the first interval cn(l) is 
decreasing, reaches a minimum and then it increases in the second interval. The minimum 
value is reached for the value of X where the spherical phonon energy is vanishing. Due 
to this behavior, in the second interval there is no breaking down for the QRPA approach. 
Due to the specific dependence of the terms 2e m and Q ^ on the strength X , cn(l) behaves 
as if from the minimum point on, the effective two body interaction changes its attractive 
nature to a repulsive one. Also, it is worth noting that in the region around the minimum 
X the difference min(2e m ) — cn(l) is large comparing it with the values corresponding to 
the A' from the complementary interval. This feature reflects the collective property 22] 
of the corresponding phonon state. We may say that the maximal collectivity is reached 
for the critical value of X. Around this point the spherical and deformed systems might be 
described in an unified fashion by using a spherical single particle basis. 
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FIG. 4: The first QRPA equation root as a function of the QQ interaction strength. The minimal 
two quasiparticle energies are also presented as a function of X. 

V. SUMMARY AND CONCLUSIONS 

In the previous sections we formulated an approach of renormalizing the QRPA such 
that no breaking down shows up. Indeed, the first QRPA energy, instead of vanishing, it 
becomes minimum and then, by increasing the long range interaction strength, is increasing. 
Things happen as if the effective interaction changes its character, from attractive to an 
repulsive one. The formalism redefines first the system ground state by accounting for the 
quasiparticle quadrupole pairing interaction. Moreover, on the top of the newly defined 
ground state a QRPA description is constructed. It turns out that the drawback of the 
standard QRPA of collapsing for a critical value of the interaction strength, is removed. In 
the new picture some higher QRPA dynamics is included. Indeed, the scattering terms are 
effectively participating in building up the new phonon operator. 

We note that the new quasiparticles are not tensors of definite rank and projection. They 
have however a definite j. This makes the difference with the picture where first one defines 
a deformed mean held and then the pairing correlations are considered. In this case j is not 
a good quantum number, but is. This difference favors the present approach, when the 
QRPA is supplemented by an angular momentum projection operation of the many body 


14 





states. 


The states considered in the present work are characterized by K = 0 and, therefore, by 
the total angular projection, finite bands of K = 0 can be defined. 
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